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Abstract. The Frobenius-Schwinger-Dyson equations are a rather high-brow abstract 
nonsense type of equations describing n-point functions of arbitrarily high composite in- 
sertions. It is not clear how to solve or even find approximate solutions of these equations 
in general, but they are worth investigating because (a certain preferred type of) renor- 
malization of composite insertions has been performed in advance: it just remains to find 
solutions given an action and renormalization conditions. Earlier work in this field involved 
only Gaussian actions or variable transformations thereof. In this work we illustrate the 
use of Frobenius-Schwinger-Dyson at a less obviously trivial level, that of the Thirring 
model. 
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1. Introduction and overview 

In a previous paper ||, we introduced the notion of Frobenius-Schwinger-Dyson equations 
associated to an action. These equations which we will remind the in next section are 
equations for renormalized expectation values (Ox, ■■, O n ) r of possibly composite insertions 
Oi. 

Under the assumption that these expectation values are obtained as 

(Ox,..,O n ) r := (R(Ox)..R(O n )), 

where (.) is a solution of the Schwinger-Dyson equation associated to the action S, and R is 
a regulator of composite insertions, and if one furthermore assumes a compatibility between 
R and S, one can derive equations for the renormalized expectation values. 

This work was financially supported by "Samenwerkingsverband FOM/SWON Mathematische Fysica", 
under grant no. MF92.3. 
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To get an idea of what this compatibility is, note the following: In case S is quadratic, 
R is typically to be though of as the normal ordering Ns associated to S, but if S is not 
quadratic then R is not to be thought of as normal ordering Nq of some quadratic action 
So, but rather as an operation which is directly related to S'.Now it is not clear how to 
generalize normal ordering to arbitrary non-quadratic actions S, so that one is forced to 
think about what it means to say that R is a regulator of composite insertions directly 
related to S. The mere requirement that infinite-dimensional calculations give finite results 
would reduce the theory to triviality since we might in that case just as well take R = 0. 
In || it was noted that normal ordering Ns for quadratic actions (or in fact a larger class 
for which it is clear how to generalize normal ordering) satisfies a compatibility condition 
with S such that 

1. The compatibility condition can be generalized to arbitrary action S. (This is fortu- 
nate) . 

2. The compatibility condition does not determine Ns once S is given. (Unfortunate). 

Forgetting about point two for a moment, we were then led to define regulators as those 
operations R that satisfy the compatibility condition, which in turn implies identities for the 
corresponding renormalized expectation values, which we called the Frobenius-Schwinger- 
Dyson equations. These equations seemed so natural that it was more or less suggested 
that computing functional integrals with composite insertions and fixed action is nothing 
but solving the FSD equations, although that point was only illustrated with Gaussian 
integrals and variable transformations thereof. 

Finally, coming back to point two one notes that Frobenius-Schwinger-Dyson equations 
associated to an action can have multiple solutions if they are not supplemented with extra 
renormalization conditions. 

In this article we will push the illustration of the FSD equations a little further by using 
them to derive the expression for renormalized 2n-point functions 

(^(x 1 )^ A Hx 2 ),^(x 3 ),...^ B ^ yi ),^(y 2 )^ B ^y 3 ),...) r 

first found by Johnson ^ for a two-dimensional Fermion field with the Thirring action. To 
fix them, we will have to impose some renormalization conditions, which we will choose 
in such a way as to remain as close as possible to Johnson's expressions. What Johnson's 
method amounts to in our setting is to guess the form of renormalization conditions that 
might allow for solutions of the FSD equations, and subsequently fixing free parameters in 
this guess just by the very imposition of these FSD equations. 

2. Reminder on the Frobenius-Schwinger-Dyson equation 
We recall the following from ||: 

Definition 2.1. Let A be a symmetric associative algebra with unit and typical elements 
f,g, let L := Der(A) with typical elements X,Y and corresponding operations fg, fX, 
Xf and [X, Y] . By a renormalized structure on A we mean the datum of S S A, together 
with extra multiplications f. r g, f- r X, X. r f and [X, r Y], also satisfying associativity, Jacobi 
identity etc., such that 

1. The . r -operations induce L = Der(A, . r ), i.e. the derivations of the multiplication 
(/><?) l— *■ f-rd <ire exactly given by the operations f t— > X. r f as X £ L. 

2. Both structures have the same unit: If = f = l. r f. 
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3. The two algebraic structures and the action S are compatible in the sense that: 

X. r (YS) - Y. r (XS) = [X, r Y]S, 

(f. r X)S = f. r (XS) - X. r f. 

The algebra A together with the renormalized structure will be called a renormalized volume 
manifold. For symmetric linear maps (, , , ) r : A® n — > K we will be interested in the following 
properties: 

1. (., .) r : A® 2 — > R is a positive non degenerate form. (Positivity). 

2. (0 1 ,0 2 ,...,O n ) r = {O x . r 2 ,..,O n ) r , (Frobenius). Q 

3. (XS,0 1 ,...,O n ) r = Y: 1 ? =1 {0 1 ,..,X. r O i ,..,O n } r . (Schwmger-Dyson). 

By substituting S/h for S and taking the limit h — > 0, we see that at H = 0, we may take 
both algebraic structures to be equal. In what follows, we will take the usual multiplication 
to have priority over the renormalized one, i.e. f '. r gh means f. r (gh). 

In this article we will concern ourselves with the case where A is the algebra of words in 
odd symbols tp A (x),^ A (x) and their derivatives, and S is given by the Thirring action. 



3. Reminder on chiral currents 

3.1.1. Chiral currents. By a chiral current we shall understand a current J which is diver- 
genceless and rotationfree if the Euler-Lagrange equations are satisfied^. More precisely, 
currents such that there are functional vector fields Y and Y au such that 



= YS, 



In the context of renormalized volume manifolds, one may now define [J^ > r .] : A — > A 
by the following properties: 



cyj^ > r O] :=Y. r O, 
d^Ju >r O] - d v [J„ > r O] := Y^. r O. 
[J^(oo),> r O] := 0. 

Note that this differential equation for [J^ > r .] can be solved rather explicitly in terms 
of an integral expression of the right hand side. 

We now see that if (J ll (oo),0) r = then the chiral currents satisfy the identity 

{J",0) r = ([J" > r 0}) r . 

Such an identity may in turn be used to produce differential equations for n-point func- 
tions in case certain identities are known for the . r operations. Such identities may be 
known in the following cases: 

1. Everything was already solved by some other method, say because S was Gaussian. 

2. One guesses properties of . r . 

A Frobenius algebra is a symmetric associative algebra with metric such that (a, be) = (ab,c). 
2 In two dimensions this reduces to the more common meaning of a current J = *J tl dx tl = fdz + gdz 
being chiral iff / is holomorphic and g antiholomorphic. 
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3.1.2. Example 1. As a first example consider the action S(<p) = J d^cf)d fl (j)d D x/2 for scalar 
fields <p : W D — > M., and renormalization condition 

o x . r — — o 2 = Oi — —-. r o 2 . 

5<j)(x) 5cj)(x) 

This renormalization condition was discussed in more detail in ||: The FSD equations for 
non degenerate Gaussian integrals subjected to this condition allow for a unique solution 
usually referred to as the solution for this action. The solution of the . r multiplication in 
terms of the usual one can for this renormalization condition and action S be inductively 
expressed as follows (see appendix, proof |l]): 

-Ol).r02 = jrniOl-rOz) + X . r -r—0 2 . 



5<j){x) 5(p(x) 8(p(x) 

Now the action above also happens to have a chiral current := d^cf), by taking Y x := 
—5/5(p(x), and := 0. The expression V p (x) := exp(pcf)(x)) then has the property that 
its derivative can be algebraically expressed in terms of j M and V p itself: 

The properties of . r that are needed to derive a differential equation for n-point functions 
of the Vp's are 

1. Y\ r Vp(y) = Y*V p (y), 

2- j^(x)V p {x) = lim y ^ x j^y). r Vp{x) - [j^y) > r V p {x)]. 

Indeed, the first equation above follows directly from the main renormalization condition; 
For the second equation we use j^(x)Vp(x) = \\m. t ^j tAi {y)V p (x) , and use the Gaussian 
recursion relation for the . r product. 

The way these two equations imply a differential equation for n-point functions is as 
follows: 

d x ,(Vp{x),0) r = p(j fl (x)V p (x),0) r =plim((j fl (y). r Vp(x)-[j ll (y) > r V p (x))),0) r 

y >x 

= p lim(j M (y), Vp(x). r O) r - > r V p (x)\, (D) r = lmi{V p {x), [j M (y) > r 0]) r 

y >x y >x 

= (V p (x), [j M (x) > r 0}) r 

3.1.3. Example 2: Next consider a slightly more complicated case, that of Fermionic spinors 
ip and "ip, with Dirac action S = J 'tpi^y^d^— im)ipd D x, together with the Fermionic analogue 
of the renormalization condition of the previous example. This case is just as solvable as 
the previous one since it is also Gaussian. The current jn := ip^^ip is chiral in case we have 
m = and D = 2: Indeed defining the functional vector fields 

8 .a, * 8 



5ip A {x) v '5^ A (x) n 

1 s s 

we have 

d^f = YS, 

d^jv - dvjv - Y^S = -{datpT^ij} - i)T pwdail) + 2z'm'i/j[7 iU ,7v]^}, 
where := 7^7,, - 7^7a7 M > 
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which is zero in D = 2. Assuming that m = and D = 2, one may again derive properties 
of . r that imply a differential equation, for re-point functions ip and ?/>'s this time: Setting 
£/j,uY '■= Y[iu, where dx 11 A dx v =: e^dx 1 A dx 2 , and using the notion 

1 3 

lim e ^ {0 0)eR 2 /(e) :=lim - ^ f(Rot ix90 o(e)), 

i=0 

so that in particular lim e _>n e a f(\e\ 2 ) = 0, and lim e _> Cpr = ^2~> we have: 

1. Y. r ip = Yip ; Y. r ip = Ytp ; Y. r tp = Yip ; Y. r tp = ftp 

2. Iim" e ^ o {.7 M 0r + e). r V' yl (x) - [j^x + e) > r ip A {x)] - j^x + e)^*(a:)} 

= ^{l u l,) A B d^ B {x). 

(See appendix, proof ||). 

3 - ^).,V^(y) = ^-y)- 

4. ^(a;) = lim^o ^ A (x). r (^ tl ) AB ^p B (x + e). 

That these properties by themselves imply a differential equation for the 2n-point func- 
tions will be shown in more generality in the next section, from which the present case is 
recovered by setting A := 0. 

4. The Thirring model with Johnson's renormalization conditions 

In the associative algebra on odd symbols %p A (x) and ip A (x) where x G M. 2 and A is an 
index in a 2D representation space for the Clifford algebra of M? with 7^ = 7^4, Thirring's 
Lagrangian is defined as ip) := tpj^d^tp + ^j^j^, where j 1 * := -ip^^ip. Consequently, we 
have 

L #r = (7 m )a^ B + Afr^AB^*. 

2. $r = fy^VW - A(7^)^i^ B . 

3. 0^ = 

4. <9^ji/ - = Y^S. 

I.e. although the Lagrangian is not quadratic any more, the current jn remains chiral even 
when A 7^ 0, which basically makes a solution of the Thirring model possible. 

Next, we need renormalization conditions to fix the renormalized 2re-point functions. 
Now it is not evident what kind of conditions can be imposed such that solutions of the 
Frobenius-Schwinger-Dyson equations satisfying them actually exist. This general problem 
will remain unaddressed in this paper: We will restrict ourselves to deriving a number of 
Johnson's results Q by using the Frobenius-Schwinger-Dyson equations: In our language, 
Johnson made a smart guess about the form of a number of renormalization conditions, 
by introducing some undetermined parameters a, a in the conditions as they hold in the 
Gaussian A = case: 

1. Y. r tp = aYip ; Y. r ip = aYxp ; Y. r tp = dYtp ; Y. r ip = aYtp, see [[[], formulae 14,15]. 

2. lhn e -> {jn( x + e). r ip A (x) - [j^x + e) > r tJj a (x)] - j^x + e)ip A (x)} 

3 We will not address the question whether this condition can be consistently extended to be a distribution 
for all values of x and y including x — y. The price we pay is that solutions of the differential equations that 
follow need not be unique. 
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4. ji*(x) = Iim e ^o Ux^a{e 2 )i) A {x). r {-i^) AB '4) B {x + e), see formula 29] and @, formula 
4.71], 

where t/^a is denned by (^ A (e), ip B {0)) X ' a > d = U^ l a ~ a (e 2 ){4) A {e) : V> B (0))^ =0 ' a=1 ' a=1 , which 
as we will see exists. We see that the Gaussian properties with canonical renormalization 
condition are recovered by setting a = a = b = 1. The parameter b does not occur in 
Johnson's work, but we will need condition |2| to make everything work. It will not enter 
the final answer because the 6-term will be multiplied by 7^, and since in two dimensions 
7 m 7^7m = 0. 

Repeating the procedure that we have seen before to produce differential equations, we 
see that the defining differential equation for [ju > r .] acting on ift, ip is solved by 

[f(x) > r ^ A (y)] = (-af»(x - y)5 A + ae^f u (x - y)i A B )i> B {y), 

[f{x) > r $ A (y)] = (o/"(s - y)S A - ~ae^f v (x - y)j A B )i> B (y), 
where 7 := 7 1 7 2 , f(x) := (In |x|)/27r, and ffi(x) := d^f(x). For example: 

d,[f(x)wi, A (y)] = -ad^{x-y)5U B {y) = -aS(x-y)S A ^ B (y) = aY^ A (y) = Y x . r ^ A (y). 

The differential equation that now holds is the following if O is a repeated . r -product of 

^(yi),^(yi)'s where x / yf. 

(hn A Bd^ B (x),0) r = -\( 7 n A B(*P B (x),[j,(x) > r 0]) r . 

Proof 



SC 

"5-i[) A (x) 



b 



= -A(7 M ) A B lim^oO; t (x + e). r ^ B (x) - {j^x + e) > r ^ B (x)} + —(Y^) B C d^ C (x), O) 
= -\{r) A B{^ B {x)i > r O}) - ^Yl,) A B(d^ C ,0). 

□ 

Up to now, renormalization condition || has not been used. We will prove in the next 
section that it implies the relations first found by Johnson 

1 . 1 



1-A/2tt 5 1 + A/2tt' 
5. Free field realization, restriction of the parameters a, a. 

The differential equation satisfied by the expectation values of the Thirring model can 
be solved by expectation values of special integrands of a Gaussian weight: Indeed, define 
the following Gaussian contractions || for two Bosonic fields </>i and 4>2 an d Dirac fields % 
and x- 

[4>i(x) 4>j(y)} := -5ijf(x - y), 

bc>x] ■= [x>x] ■= 0, 
[x A (x)>x B (y)\ ■■= dn AB U(x-y) =■■ [x A (x)> x B (y)}. 

Following [Q, formula IV. 5], let us define 

^ A (x) := ( e +XklMx)+Xk2Mx ^) A B x B (x), 
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for yet to be fixed numbers fei, k 2 - Making the Dirac operator j^d^ act on these expressions 
will bring down the matrix 7, and using that 7^7 = e^ju, we become interested in 

f ± := hd^x ± k 2 e^d v <h 

By Gaussian Wick calculus we then have: 

[£(x) > * A (y)] = {Xklr(x - y)5 A - \k\e^ f u (x - y)j A B }* B (y), 

[f_(x) > V A (y)] = {-Xklr(x - y)5 A + \k 2 2 e^f u (x - y)7 A B }* B (y), 

= => h,) A B (d^ B (x), 0) r = -\h») A B(* B (x), [£(*) > 0]) r . 

The last differential equation is exactly the one derived using Johnson's renormalization 
conditions, provided the action [jt > .] equals the action of j 7 ', which happens if we have 
— Xk\ = a and —\k\ = a. Consequently the n-point functions of the ^'s above with these 
values of fcj satisfy the differential equation for the ^'s of the Thirring re-point functions. 
In particular the two-point function can be computed, leading to 

(i; A (e), 4> B (0))$' a ' d = e A(a - a) ^ (0 ^ (e)) (V; j4 (e),^(0))^ =0 ' a=1 '' i=1 , 

so that [/^ 5 (e 2 ) = e x ^ a ~ a ^^ ^^ . Furthermore, in view of renormalization condition || 
the Thirring expectation values of j^x) are reproduced by the following expression in terms 
of the Gaussian fields: 

J"[x) := lirn^o U x ^(e 2 )^> A (x). r ( 7ll ) AB ^ B (x + e). 

But using Gaussian Wick calculus, one may prove (see appendix, proof ||) that 

= <s> A mAB* B + ^# = x^Ubx 3 + ^# =■■ f x + ^jI- 

Since 0) T = ([j^ o r 0]} r , and comparing > r ip A ] with \jg + -^j^ > \t A ], we see that 
the following two expressions must be equal: 

1. (-af(x - y)8 A + ae^Mx - y)j A B )^ E \y) , 

2. (-(1 + %)f*(x - y)5 A + (1 - §)s^U(x - y)j A B )Tp B (y), 

giving that a = 1 _1/ 2 tt ® = i+a/?tt as P rormse d- 

6. Partial vs. Total renormalization conditions. 

By partial renormalization conditions we shall understand renormalization conditions 
which together with the FSD equation do not completely determine the . r operations, as 
opposed to total renormalization conditions which will be those which do fix the . r 's: For 
example fdi. r g = f. r dig is a total renormalization condition for Gaussian actions. 

Johnson's renormalization conditions for the Thirring model seem to be only partial, since 
when restricted to the Gaussian A = case, Johnson's conditions are obviously weaker than 
fdi- r g = f-rdig. 

One may ask if there is a natural way to strengthen Johnson's renormalization conditions 
so that they also determine the operator product or expectation values of higher order 
insertions. A simple try at guessing such an extension is to first extend the free field 
realization F : ip 1— > exp(— Xk\4>i — Xk2<p2^)ipo by F(ab) := F(a)F(b), and then trying to 
pull-back the Gaussian . r operations by F. This however fails because F(a). r F(b) is not in 
the image of F for all a and b. Another try could be to use the bosonization procedure as 
a free field realization. The same problems as above might in that case occur although we 
didn't check that. 



s 
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7. Conclusion 

What we have seen is that the Probenius-Schwinger-Dyson equations can reproduce 
known expressions for renormalized expectation values, the advantage being that there 
was no need of explicitly regularizing and renormalizing the fields. It remains to be proved 
that there actually exists a solution of the FSD equations for the Thirring model. 

It would seem that the same method will apply for example to the Wess-Zumino-Witten 
model, since the derivation || of differential equations for n-point functions in that case 
is not different from Johnson's derivation: Arbitrary multiplicative parameters k = R have 
been introduced H], formula 3.1], and consistency with renormalization conditions on higher 
composite operators fixes the value of k in terms of the integer coupling constant k of the 
model H formula 3.11]. 

^From a more fundamental point of view it remains an open question, even in finite 
dimensions, whether there is an algorithm which given an action and renormalization con- 
ditions will produce approximations of solutions of the FSD equations up to arbitrarily 
high precision. Unlike approximation schemes for positive solutions of the usual Schwinger- 
Dyson equation, which are just the approximation schemes for integration, an approxima- 
tion scheme for solutions of the FSD equations seems difficult to construct given that the 
number of solutions will in general depend on the renormalization conditions. 
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Appendix A. Proofs 

1. The renormalization condition f. r dig = fdi. r g implies (f. r g)di = f- r (gdi) and {{diS)0\). r 02 
(a j 5)(O 1 . r O 2 ) + O 1 . r i O 2 : 

Proof 

To prove the first formula, we have to show that both sides of the equation have the 
same . r action: 

[(f. r g)d]. r h = (f. r g). r dh = f. r {g. r dh) = f. r (gd. r h) = (f. r gd). r h. 

The second formula now follows from the first as follows: In the formula f- r XS — 
X. T f = (f-rX)S that holds by definition of renormalized volume manifolds, take 
/ := e> 2 , and X := O x di- This gives 

2 . r (Oid t S) - O x . r di0 2 = Qi-riOxdiS) - O x di. r 2 

= f. r XS - X. r f = (f. r X)S = ^.rO^S 

By the first formula, the last term equals (C^-rOi^S • 

□ 

2. For Gaussian Dirac fields we have: 

hm e _^o{i M (^ + e). r if; A (x) - [j^x + e) > r i/) A (x)} - j^x + e)i> A {x)} 

= ^(i v i,) A B d^ B {x). 

Proof 

Indeed, the left hand side is the limit of 

[iP A {x) o j^x + e)] - [j^x + e) > r ^ A (x)} 
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= [^(x) >i, B (x + e)]( 7M W C (* + e) + f v (e)(l v 1 ») Aci>° (x) 
= -(Y^)Acf P (e){^ C (x + e) - V C (x)} = -{Yl»)Ac^d^ c {x) + 0(e), 
which in the limit approaches right hand side. 

□ 

3. The different currents introduced in section [| satisfy the following identity: 

j» = * A (jnAB* B + = x A (r)ABx B + ^# =: % + 

Proof 

First by expanding roughly speaking ^ = e~("^x an d ^ = e + ^'"''x> so that ^f. r ^f = 
(e + x)-r(e _ x) = {e + }. r {e"}(xx + (xx)); an d using the matrix version of e px . r e qx = 

e pq e (p+q)x^ we have . 

e -Ao(0(a:)0(j/))-Aa(0(a!)0(i/))7*7»^/ a .N r ^r(y)) 

where 7 X acts on i- e - (7xX( a; )x(y)) A ' B = (7 j4 CX C "( 3; )x S (?/))• «LFrom this formula 

one deduces the identity among currents by taking the limit y — > x: 

#(x) = ^ A (x) 7 ^^ B (x) = linT^o * A (x)^ AB y B (x + e=:y)=I-II, 

where with B e := ((f)(0) (j)(e)): 

I = m e ^(e~ XaB ^ x ~ aB ^y) A c B D^ c (x). r ^ D (y)i% 

= m e ^ e- x ^ B ^ c (x). rl ^ D (y) = J*(x) 
II = lim" e _, [ e ^^(^( i O-0lfe))+v^(^(x)7,-02(^/)7 H )]A c B D ( 7 P)C7D^^ _ y )^) AB 

= IEi e ^o e V ^ {4>l{x) ~ <l,l{y)) fp(x - y )Tr(e' / ^^' 2( - x) ~^ ( - y) ^Yl fl ) 

x{2 5 ^cos[v /Z A^(0i(x) - Mv))\ ~ 2e"*sia[V=A^i(:r) - MvW 
ImT^o [1 - V^e^^x) + 0( e 2 )]_^2[^(l + 0(e 2 )) + £^v^e CT c^ 2 (x) + 0(e 2 )] 

= ^{V=ted P Mx)g pl1 - e^V^xbldpMx)) = ^f + (x) 

□ 
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